clear all;close all; clc; 
r=1.25;  % сопротивление для диода 2Д103А из примера
C=33e-12; % емкость конденсатора задана по условию
C1=C/2; % последовательное соединение двух конденсаторов одинаковой емкости
tau=r*C1;
fprintf('tau=%f ns\n', tau*1e9); % вместо %f записывается то, что стоит после запятой
T=0.01*tau; % Шаг
t=0:T:3*tau; % Создать массив с шагом T от 0 до 3 tau
lt=length(t); 

% выделение памяти

Uc=nan(1,lt);


Uc1(1)=0;    % начальное значение выходного напряжения на конденсаторе
            % при отсутствии входного воздействия


 % Воздействие на скачок
S=2;  %постоянное входное воздействие равное 2
alpha=tau/T;
for k=2:length(t);
    Uc1(k)=(abs(S)+alpha*Uc1(k-1))/(1+alpha);
    Uc(k)=Uc1(k)/2;
    end;
    figure(1)
plot(t,Uc);

%   Построение АЧХ



% Воздействие в виде двух функций Хевисайда
for k=2:lt;
    if k<lt/2;
        S(k)=1;
    else S(k)=0;
    end
    Uc1(k)=(abs(S(k))+alpha*Uc1(k-1))/(1+alpha);  
    Uc(k)=Uc1(k)/2;   %  выходной отклик
    end;
    figure(2)
plot(t,Uc);

f=0:1e6:1e12;
omega=2*pi*f;
K=1./(sqrt(1+(omega*C1*r).^2));  %  коэффициент передачи RC-цепи
figure(3)
plot(omega,K);



